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^3 : I. INTRODUCTION 

The anisotropic harmonic oscillator potential is of physical interest in quantum mechanics, 
^ ' since it may describe the motion of, say, an electron in an anisotropic metal lattice. With 
■ external electromagnetic fields imposed on, such models may also be useful in semiconductor 
physics [1,2]. For an isotropic harmonic oscillator potential and a homogeneous static mag- 
netic field the Schrodinger equation can be easily solved in the cylindrical coordinates. If 
the harmonic oscillator potential is anisotropic, the problem is not so easy. Nevertheless, it 
can be solved analytically because the Hamiltonian is quadratic in the canonical coordinates 
and momenta. There exists quite general formalisms for solving systems with such quadratic 
Hamiltonians [3], and the above problem represents a typical example [4]. In Ref. [4] only the 
isotropic case is considered as an example for illustrating the general method. The anisotropic 
case with a homogeneous static magnetic field was studied in Refs. [5,6]. In these works the 
magnetic field was arranged to point in one of the rectangular axis directions. In this paper we 
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will solve the problem in a somewhat easier formalism. A homogeneous static electric field is 
also included, though its efi^ect is trivial. We will also deal with a case where the magnetic field 
has a more general direction. This seems not having been considered before. Coherent and 
squeezed states of these systems will be discussed in some detail. In terms of these concepts 
the time evolution of some wave packets can be discussed very conveniently. 

In the next section we consider a charged harmonic oscillator with arbitrary frequencies 
ujxi ujy and lv^ in three rectangular axis directions. A homogeneous static magnetic field in 
the z direction and a homogeneous static electric field in an arbitrary direction are imposed 
to the system. The problem can be easily reduced to one equivalent to the case without 
an electric field by making some coordinate transformation. In Sec. Ill we diagonalize the 
reduced Hamiltonian and give the energy eigenvalues. Although this has been previously done, 
the method used here is different and seems more convenient, especially for the subsequent 
discussions of coherent and squeezed states. The wave functions are explicitly worked out in 
Sec. IV. In Sec. V wc discuss some special cases, some of which may need different handling, 
or even cannot be solved. In Sec. VI we consider a somewhat different case where = 
and Ez = {Ez is the z component of the electric field), but the magnetic field has a x or y 
component, in addition to the z component. To our knowledge, this was not considered before. 
It is solved analytically by the formalism of Sec. III. In Sec. VII we discuss the coherent and 
squeezed states and their time evolution in detail. A brief summary is given in Sec. VIII. 



II. REDUCTION OF THE HAMILTONIAN 

Consider a charged particle with charge q and mass M, moving in an anisotropic harmonic 
oscillator potential, a homogeneous static magnetic field B = Bsz, and a homogeneous static 
electric field E = E^e^ + EySy + E^e^, where e^,, ey and are unit vectors of the rectangular 
coordinate system, and B, E^, Ey and E^ are all constants. The stationary Schrodinger 
equation is 

//^*(x) =£:^*(x), (1) 

where the Hamiltonian is 

= ^ (P - ' a)V '-Micoy + uy + coy) - ,E . X, (2) 

where E ■ x = E^x + Eyy + E^z. Since the magnetic field points in the z direction, one may 
take Ax and Ay independent of z, and Az = 0. Then the Hamiltonian can be decomposed as 

+ ^Miuy + uy) - q{E,x + Eyy), (3) 
and 

Therefore the total wave function can be factorized as ^(x) = ip{x, y)Z{z), where y) and 
Z{z) satisfy the following equations: 
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(5) 



H'Z{z) = S'Z{z), (6) 
with -\- £^ — . Obviously, can be recast in the form 

where zq — qEz/Mul, so that Eq. (6) can be easily solved with the following results. 

^"3 = (^3 + ^) /^o;. - n3 = 0,l,2,..., (8) 

Zm{z) = - ^o), na = 0, 1, 2, . . . , (9) 

where we have quoted the standard wave functions for harmonic oscillators: 

C(^) - , exp , n = 0,l,2,.... (10) 



Therefore our main task is to solve Eq. (5). We take the gauge 

A, = -\By, Ay = \Bx, (11) 
where x — x — xq^ y — y — yo, and xq — qEx/MujI, yo — qEy/McUy, then if^^ can be written as 

- ^(PJ + P^) + ^ '-M.lf - ..i. - 15 - (12) 

X y 

where = xpy — yp^, ujb = qB/2Mc, and 

^1 = \/^x+^B, ^2 = y/ujf-H^. (13) 

Note that a;^ may be either positive or negative, depending on the signs of q and B. Now Eq. 
(5) can be written as 

H'ilj{x,y) = £ilj{x,y), (14a) 

where 

„2 7712 „2 Tp2 

^"^ +2Ma;2 + 2Ma;2- ^l^) 

Eq. (14) has the same form as that for an anisotropic harmonic oscillator moving on the xy 
plane under the influence of a homogeneous static magnetic field perpendicular to the plane, 
except that x, y are replaced by x, y. The main feature is that the reduced Hamiltonian is 
quadratic in the canonical variables. Though this equation has been studied by some authors, 
we will solve it in a somewhat different and simpler way in the next two sections. Here we 
point out that the calculations below will be simpler if one take the gauge A^ — —By, Ay — 
or Ax — 0, Ay — Bx. However, we prefer the above gauge (11), since in this gauge it would 
be convenient to compare the results with known ones when cux — ^y-i or uji — loi. 



III. DIAGONALIZATION OF THE REDUCED HAMILTONIAN 



If cux = LUy, Eq. (14) can be easily solved in the cylindrical coordinates. In the general case 
this does not work, so other methods are necessary. We define a column vector 

X ^ (x,Px,y,Pyy, (16) 

where r denotes matrix transposition, then the reduced Hamiltonian (14b) can be written as 

H = ^x^nx = ^x^nx, (17) 

where is a c-number matrix which we take to be symmetric. We do not write down Tl 
here since this is easy. If one can diagonalize Ti. then the Hamiltonian would become a sum 
of two one-dimensional Hamiltonians for harmonic oscillators or something similar (repulsive 
harmonic oscillators or free particles). However, a crucial point here is to preserve the com- 
mutation relation 



[X^,Xp]^-h{J:yU, a, 1,2,3,4: (18) 

after the linear transformation that diagonalizes Ti, where Sj^ = diag((T,y, (jj,), and ay is the 
second Pauli matrix. Therefore the needed transformations are the so called symplectic ones 
[3,4]. Here we will transform the Hamiltonian into a form expressed in terms of raising and 
lowering operators, so the formalism is somewhat different. Let us proceed as follows. 
Because H is quadratic in the canonical variables X, it is easy to show that 

[iH,x] = hnx, (19) 

where 



(20) 



is another c-number matrix. As both 7i and illy are real, so is fl. Our first step is to diagonalize 
so we write down it here. 






1/M 







-Muj 








ujb 


-ujb 








1/M 





-UJb 


-MujI 






y u — UJB — 1V1IJJ2 "J J 

The characteristic polynomial for this matrix is 

det(A/ - fi) = + + c, 

where 

b = ujl + ujl + AujI, c = ujIujI- 

Since 



[21) 



(22a) 



(22b) 



A 



A = b'-ic={u;l- ulf + Suliul +u;l + 2a;|) > 0, 



(23) 



we have two real roots for A^. If cUa; = or ujy = 0, one of the two roots is zero, and the 
following discussions are not valid. We will return to this case latter. Currently we assume 
that both ujx and ujy are nonzero, then c > and -\/A < 6, and both roots for are negative. 
Therefore the above characteristic polynomial has four pure imaginary roots: 

{Ai, A2, A3, A4} = {-icTi, icTi, -ia2, ia2}, (24) 

where 

and (Ti > o"2 > 0. The equal sign appears when ujb = and = Uy. 

Because fl is not symmetric, right eigenvectors and left ones are different. Wc define two 
left eigenvectors (row vectors) Ui corresponding to the eigenvalues — icTj (i = 1, 2) by 

UiD, = —iaiUi, i = 1, 2, (26) 

then the other two are u*, corresponding to the eigenvalues iaj {i = 1, 2). Similarly, the right 
eigenvectors (column vectors) are Vi and v* {i = 1, 2), satisfying the equations 

flVi = -iaiVi, i = l,2 (27) 

and the complex conjugate ones. Prom these eigenvalue equations it can be shown that 

u;vj = UiV* = 0, yt,j = l,2, (28a) 

and UiVj = u*v* = for (jj 7^ aj (this should be checked individually when ai = (72, see below). 
By appropriately choosing the normalization constants we have 

UiVj = u*v* = 6ij, W i,j = 1, 2. (28b) 

The specific expression for the above eigenvectors are not necessary in this section. However, 
a relation between the left and right eigenvectors is very useful in the following. Using Eqs. 
(20), (26) and (27), it is easy to show that Vi oc so we choose 

Vi = -^yul i = 1, 2. (29) 

Note that the condition (28b) only determine the product of the normalization constants in 
Ui and Vi. With this relation both constants are fixed (up to a phase factor). Also note that 
Eqs. (28b) and (29) lead to —UiT^yuj = 1 (no summation over i on the left-hand side). The 
left-hand side of this equation is real (which can be easily shown), but not necessarily positive. 
Actually it can be shown by using the eigenvalue equation that —UiHyu} = {ui'Ey)li.{uiT^y)y ai. 
Thus the sign of this quantity depends on the matrix 7i, and in general Eq. (29) should be 
replaced by t'j = ejSyM| where ej = ±1. For the present case, however, Eq. (29) is sound. 
Now we define a 4 x 4 matrix Q by arranging the column vectors in the following order: 

Q = {vi,vl,V2,vl). (30) 



Using Eq. (28) it is easy to show that 

g-i = K,<,«5,<r- (31) 

With the help of the eigenvalue equations (26) and (27), we have 

Q~^QQ = diag(-icri, ifTi, -ia2, i(J2). (32) 

Therefore the matrix fl is diagonalized. Using the relation (29), it is not difficult to show that 

= -^,Q~%, (33) 

where = diag((T2, cr^), and is the third Pauli matrix. This relation is very useful in the 
following. 

Next we will diagonalize the Hamiltonian. With the above results this is easy. We define 
two operators 

Oj = UiX/ y/h, i = 1,2. (34a) 

Their hermitian conjugates are 

a| = u*X/Vh, i = 1, 2. (34b) 

Using Eqs. (18), (28) and (29) it is easy to show that the nonvanishing commutators among 
these operators are 

[ai,a]]^5ij, i,j = l,2. (35) 

Therefore they are similar to the raising and lowering operators for harmonic oscillators. We 
define a column vector 

{ai,a\,a2,aly, (36) 
then Eq. (34) can be written in the matrix form 

A = Q-'X/Vh. (34') 

The inverse is 

X = VhQA, (37) 

and thus 

= X^ = VUA^Ql (38) 

Substituting these into Eq. (17), and using Eqs. (20), (32) and (33), it is rather easy to show 
that 

H = |MtE^, (39) 
where S = diag(cri, cti, cr2, £72). Using the commutation relations (35) this becomes 

H = hai{a\ai + |) + ha2{a\a2 + |). (40) 



Therefore the Hamiltonian is diagonahzed. More specifically, it becomes the sum of two one- 
dimensional harmonic oscillator Hamiltonians. The energy levels are readily available: 

^nma = + +^cr2(n2 + I), ni,n2 = 0, 1,2, . . . . (41) 

When ojb = 0, we have cti = max(c<Ja,, ujy) and (T2 = vcnxi^uJx, ojy), and these energy levels reduce 
to known results, as expected. Substituting these into Eq. (15) we obtain the energy levels for 
the motion on the xy plane: 

4'f„ = (m + i) + n.. (". + 5) - S - 0. 1. 2, . . . . (42) 

The total energy is the sum of £nxn2 ^"^^ ^nz given by Eq. (8). The wave functions on the xy 
plane will be worked out in the next section. 

The method employed in this section is similar to that used in Ref. [7] which dealt with a 
charged particle moving in a rotating magnetic field (also studied in Refs. [8] and [9]). However, 
we present here a relation between the right and left eigenvectors, which leads to the very useful 
relation (33), so that the calculations here seems more straightforward and easier. The method 
can be easily extended to systems with more degrees of freedom. 



IV. THE WAVE FUNCTIONS IN THE COORDINATE REPRESENTATION 

The eigenstates of the Hamiltonian (40) corresponding to the eigenvalues (41) are 

1^1712) = ^i==(al)"H4)"'|00), 711,712 = 0,1,2,..., (43) 

where |00) is the ground state satisfying 

ai|00) = aslOO) = 0. (44) 

The task of this section is to work out these states in the coordinate representation, that is, 
the wave functions 

i>mn2{x,y) = {xy\nin2). (45) 
For this purpose the left eigenvectors are necessary. These are given by 

Ui = Ki{ - iMai{a^ - - 2a;|), - Musia^ + o;^), i2a;sa0, i = 1, 2 (46a) 
and their complex conjugates, where the normalization constants are 

= {2Ma,[ia^ - + ^^W^Y'I\ (46b) 

When ijJb = 0, {(Ti\i = 1,2} = {uJx,ujy}. For Ux, the last two components of Ui vanish. For 
LVy, the first two components of Ui vanish. Thus the orthogonality is ensured. If the further 
condition uj^ = holds (which yields cxi = (J2), one may take the limit ujy — > uj^ after 
everything is worked out. 

First we should work out the ground state. In the coordinate representation, Eq. (44) 
takes the form 
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{^ijXj - ihr)ijdj)ipQo{x, y) = 0, i = 1, 2, (47) 
where Xi — x, X2 — y and similarly for xi and X2, dj — d/dxj — d/dxj, and we have defined 



two 2x2 matrices 



Mil t<13 



V 



Ui2 Uu 



U22 U2A 



(48) 



where Wj/j (i = 1, 2; /3 = 1, 2, 3, 4) is the /3th component of Uj. Suppose that 

"000(2;, y) = A^oexp[-s(x,y)], s(x,y) = Xik^jX^l^h, (49) 

where ^S" is a 2 x 2 symmetric matrix whose elements are complex numbers, and A^o is a 
normalization constant. Substituting this into Eq. (47) we obtain 

A = \^-\. (50) 

This can be easily worked out. The elements are 



All 



MLo^{ax^a 2) _ ^^2 



A 



Mujy{ai + £72) 



22 



y 



A 



12 



A21 — = lAlA 



UJx +00: 



(51a) 
(51b) 
(51c) 



y 



where the relation aia2 = uj^^Jy has been used. Therefore the ground-state wave function is 



^ooix,y) = Noexp{-lXlx^ - ^Xly^ - iXxyXy), 
where the normalization constant is 

1/2 



XxXy 

TT 



M{ai + 0-2)^/07; 



(52a) 



(52b) 



7rn(LJx + UJy) 

When cub — 0, we have = Mujx/Ti, Xy — ^^Mcuy/h, and Xxy — 0. On the other hand, if 

UJx — ujy, we have uji — UJ2, ai — u!i + \u!b\, (T2—0J1 — \ojb\-i and A^; = = Muji/h, Xxy — 0. 
These are all expected results. 

It is remarkable that the uncertainty relation for the above ground state is 



\2 ^ 



1/2 



h 

> -. 

- 2 



(53) 



The equal sign holds only when B = or ujx = ujy. Therefore the ground state is in general not 
a minimum wave packet, though it is of the Gaussian type. It seems that this is not noticed 
in the previous literature. 

For any function F{x, y), it is easy to show that 

i^ijXj - ihri*jdj)F{x, y) = -i;iexp[s*(a;, y)]r]*jdj{cxp[-s*{x, y)]F{x, y)}, i = 1, 2. (54) 
According to Eq. (43), and using the above relation we obtain 

i;r^,n,{^, y) = ^ ^J-L_^ N^K^' exp^A^x^ + iA^y^ _ iA^^^y) 

[{al - u^dx - i2ujBa,dyr[{al - u^dx - i2ujBa2dyr exp{-Xlx' - X^). (55) 



s 



V. SOME SPECIAL CASES 



In this section we discuss some special cases where some of the parameters are zero. 

1. If a;^ = and £'2 = 0, the motion in the z direction is free, and > is continuous. 

2. If a;^ = but 7^ 0, the motion in the z direction is that in a homogeneous electric 
field, which is also well known, and is also continuous, and can take on any real value [10]. 

3. In the preceding sections we have assumed that both uj^ and LOy are nonzero. In this 
case, if any one of E^ or Ey is zero, it can be substituted into the above results directly and 
no modification to the formalism is needed. However, if one of lo^ and ujy is zero, the situation 
is different. We would deal with a;^ = in the following. The other case could be discussed in 
a similar way. 

4. \i LOy — ^ but Ey 7^ 0, the Hamiltonian could not be reduced to a quadratic form, and 
to our knowledge the problem could not be solved analytically. 

5. VLijjy — ^ and Ey = 0, we choose the gauge = 0, ^4^^ = Bx, and let 



i^{x,y) = ^=e^])(iky)(p{x), (56) 



then (f){x) satisfies the equation 



^pl + \Mil.l{x-x,f^ 



(f){x) = <S^0(x), (57) 



where 



^1 = \J^l + = , (58) 



and 

{qE,^ + 2khujBf h'^k^ 
2MCjI ~2M 

The energy levels and wave functions for Eq. (57) are readily available, so the final results are 

{qE^ + 2khujBY ^ 
2J 2MCjI ^ 2M 

and the corresponding wave functions are 



' oS/r,~,2 os/f J 



O = n^An,+ '-]- > + n,=0,l,2,...,^G(-oo,+oo), (60) 



V'niifc(2;, y) = e-x.^{\ky)'il^':^\{x - Xk). (61) 

V ^TT 



If further uj^. or E^ or both are zero, the results can be obtained from the above ones directly. 

6. Finally, if both uj^ and ujy are zero, one can appropriately choose the coordinate axes of 
the xy plane such that Ey — 0. Thus it is a special case of the above case 5. 



VI. A MORE GENERAL MAGNETIC FIELD 

In this section we set cu^ — and E^ — 0, and consider a magnetic field with an x or y 
component (but not both) in addition to the z component. To our knowledge this was not 
considered previously. Without loss of generality we take 



o 



B = E'e^ + Be^. (62) 
The gauge convenient for the present case might be 

A, = 0, Ay = Bx = B{x-Xo), A, = B'y, (63) 
where Xq is the same as defined before. In this gauge the total Hamiltonian takes the form 

11 11 B? 

" m^^'^^y^ + 2M^^' ~ 2Ma;^y)2 + -Mu^x^ + -M^y - qE,y - 2uBxpy - 

(64) 

where ujb and cDi are the same as defined before, and uj'q = qB'/2Mc. Pz is obviously a 
conserved quantity, so we make the factorization 

^(x) = -j= exjp{ikz)i/j{x,y). (65) 



The Schrodinger equation (1) is then reduced to the form (14a) where now 

H-^ipl+ pD + \m^Ix^ + \M^lf - 2ujBxpy, (66) 



and 



q^El {qEy + 2fc^a;^)^ _ 
2Muol ^ 2MCol ~2M 



In these equations y — y — Vk which is different from the previous one, and 

qEy + 2kUuj'Q 



u^^^u^l + Au'l yk="-^^j^^^- (68) 

Now the reduced Hamiltonian H is quadratic, and different from Eq. (14b) only in the last 
term. Thus the problem can be solved in much the same way as before. We only give the 
results here. The energy levels are 

^ ( ^\ ^ f 1\ (l^El {qEy + 2knu'Bf h^k^ 

77,1, n2 = 0, 1, 2, . . . , /c e (— oo, +oo), (69) 

where cxi and a2 are given by Eqs. (25), (22b) and (23), but with ujy replaced by uj2- The wave 
functions are given by 

1 

*nin2fc(x) = -j=exp{ikz)ipmn2kix,y), (70) 
V 27r 

where ipnin2k{x, y) has the same form as Eq. (55), with ujy replaced by 0)2 everywhere, including 
in the expressions for A^;, \y, Ki and K2, but here 

_ 2Mu;Boj2 
n[ujx + UJ2) 

Also note that y depends on k, which is the reason why ipnin2k{x,y) has the subscript k. The 
orthonormal relation is 

/ *n;nife'(x)*nin2fe(x) dx = 5{k - k')5n,n',Sn2n',. (72) 

By the way, the wave functions (61) satisfy a similar orthonormal relation. 

As before, the ground state obtained here, though being of the Gaussian type, is not a 
minimum wave packet, except when B — (then B = E'e^)- 

1 n 



VII. COHERENT AND SQUEEZED STATES 



Coherent and squeezed states are useful objects in quantum mechanics and quantum optics. 
These are widely discussed in the literature. In the presence of magnetic fields, some discussions 
can be found in Ref. [11]. There are also examples on other apphcations of such states [12]. 
However, such states for anisotropic harmonic oscillators in the presence of magnetic fields 
are not discussed previously, to the best of our knowledge. In this section we will discuss the 
definition, the properties and the time evolution of such states for the case treated in Sec. Ill 
and IV. The discussions can be easily extended to the case of Sec. VI. 

Since the motion in the z direction is separable and is essentially that of a usual harmonic 
oscillator, we only discuss the coherent and squeezed states on the xy plane (which is a simple 
translation of the xy plane) , and their time evolution as governed by the Hamiltonian H given 
by Eq. (14a) or (40) (the additional constants in H^'^ has only trivial consequence for the 
time evolution). The definitions used below are natural generalizations of those for a single 
harmonic oscillator [13-15]. 

We define a unitary displacement operator D{pL\^0L2) = -Di(ai)D2(«2) where 

Di{ai) = ex.-p{aial - a*ai), i = 1, 2 (73) 

and the aiS are complex numbers. For an arbitrary state \ip) one can define a corresponding 
displaced state 

\(p,aia2)D = D{ai,a2)\(p) = Di{ai)D2{a2)\(p). (74) 
It is easy to show that 

D\ai, a2)aiD{ai, a2) = ai + ai, i = 1,2, (75) 

so if |</7) = 1 00) is the ground state we have 

ai|00, aia2)D = ttilOO, aia2)D, i = 1,2. (76) 

Thus the coherent states may be defined as the displaced ground state |00, q;iq;2)i3. In terms 
of the original variables the displacement operator takes the form 

D{ai, a2) = exp (^-^^f ^xp exp (^-^^FPij > (^7) 

where 

xf = iVh{ajr]*i-c.c.), pf = -iVhiaj^i - c.c). (78) 

If the wave function for the state \(p) is (f{x, y), then the one for the displaced state |</7, aia2) d 
is 

Thus we see the reason why the operator D{ai,a2) is called a displacement operator. Prom 
this and Eq. (52) it is easy to obtain the wave function for the coherent states. Obviously the 
displaced state and the original one have the same shape in the configuration space, thus the 

1 1 



uncertainty AXa in |v?, aiQ;2)i3 is the same as that in \(fi). This can also be easily shown by 
using Eqs. (37) and (75). 

Now consider the time evolution of the displaced states. It is easy to show that 

e-'"^l^aie'^*/^ = exp(i(7ii)ai, i = 1, 2. (80) 

If the state at the initial time t = is \4>{S^)) = q;iQ;2)d, then the state at the time t is 

|V(i)) = \<p{t),aua2t)D = D^{au)D2{a2t)\<p{t)), (81) 

where = exp(— io".jt)ai and \^{t)) = e^'^*'''*|(/9). Therefore if \^{f)) is known, |V'(^)) can be 
obtained by a time-dependent displacement. A simple special case is \^) = \n1n2), that is 

\tp{0)) ^ \nin2,aia2)D, (82) 

a displaced number state. In this case \(p{t)) = exp[— i(ni + |)crit — i(n2 + ^)o'2t]\nin2), and 

\ip{t)) = exp[-i(ni + l)ait - 1(^2 + |)(72t]|nin2, aita2t)D- (83) 

This means that |'0(t)) is also a displaced number state, except that the displacement pa- 
rameters are time dependent, and a time-dependent phase factor is gained. The position and 
velocity of the wave packet ■0(t,x,y) = {xy\ip{t)) is characterized by the expectation values 
xl{t) = {i/j{t)\X^\'il;{t)). Using Eq. (37) it is easy to show that 

^a(^) ~ \/h[exp{—iajt)ajVja + c.c.]. (84) 
By straightforward calculations one can show that they satisfy the following equation 

x:,{t)=Q^px'^{t), (85) 

where the eigenvalue equation (27) has been used, and at t = they give the same results as 
those obtained by using Eqs. (78) and (79). On the other hand, the equations of motion for 
Xa governed by the Hamiltonian H according to classical mechanics are 

— {Xa, H}pB — QafsXp, (86) 

where the Poisson bracket {Xa, H}pb turns out to be of the same form as the commutator 
[Xa,H]/ih because H is quadratic in X^. The above two equations mean that the center of 
the wave packet moves like a classical particle. In fact, it can be shown that this holds for 
any wave packet in any quadratic system. For the displaced number states, some more specific 
properties should be emphasized. First, the shape of the wave packet keep unchanged with 
time. Second, the center of the wave packet is oscillating, but the motion is in general not 
periodic except when a2/(Ji is a rational number. 

Next we discuss the squeezed states. We define a unitary squeeze operator S{(iX2) — 
-Si(Ci)'52(C2) where 

SiiCi) = expilQalal - K>i«i)' ^ = 1> 2, (87) 

and the Ci's are complex numbers. For an arbitrary state \(fi) the corresponding squeezed state 
may be defined as 
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!</:', ClC2)5 = 5(Ci, C2)k) =^l(Cl)^2(C2)|<^). (88) 

The following equation is useful for subsequent calculations. 

^^(Ci, C2)a^S{Cl, C2) = ai cosh Pi + aje'"^* sinhp^, i - 1, 2, (89) 

where Pi = \Ci\ and 0i = argC^. 

Unlike the displaced state, it is difficult to obtain the explicit wave function for the squeezed 
state in terms of the wave function for the original one. However, in the squeezed number state 
\n1n2, CiC2)5, it can be shown by using Eq. (89) that 

{Xa) = 0, (90) 

AXa = {h{2nj + 1) [\vja\'^ cosh 2pj + Re{v]^e"^^) sinh 2pj]y/^. (91) 

Compared with the results for — (2 — 0, we see that the center of the squeezed state is the 
same as that of the original one, but the uncertainties are different. Thus the states are indeed 
"squeezed" . 

Now consider the time evolution of the squeezed states. If the state at the initial time t = 
is IV'(O)) = Iv?, CiC2)5, then the state at the time t is 

= W)XuC2t)s = ^l(Clt)52(C2*)|^(t)), (92) 

where Qt = exp(— i2(Tjt)^j and |(/5(t)) = e~^^*^^\ip). Therefore if |v5(t)) is known, \ip{t)) can be 
obtained by a time-dependent squeeze. A simple special case is \ip) = \n1n2), that is 

1^(0)) = |nin2, (1(2)5, (93) 

a squeezed number state. In this case \(p{t)) — exp[— i(ni + |)crii — i(n2 + ^)cr2t]\nin2) , and 

\ijj{t)) = exp[-i(ni + |)crit - i(n2 + ^)(T2t\\nin2,CitC2t)s- (94) 

This means that \ip{t)) is also a squeezed number state, except that the squeeze parameters 
are time dependent, and a time-dependent phase factor is gained. Though it is difficult to 
obtain the wave function ip{t, x, y) — {xy\ip{t)) explicitly, in these states it can be shown that 

{Xa)t = 0, (95) 

and 

AtXa = {h{2nj + l)[\vjaf cosh 2pj + Re{v^^e''*'^-'^''^*) sinh 2pj]}^/^ (96) 

These results mean that the center of the wave packet is at rest, but the uncertainties are 
oscillating. Thus the shape of the wave packet changes with time apparently. As before, the 
motion is in general not periodic except when a2/<Ji is a rational number. 



VIII. SUMMARY 

In this paper we have studied a charged anisotropic harmonic oscillator moving in a ho- 
mogeneous static electromagnetic field. Several configurations of the electromagnetic field are 
considered. One of these configurations has been studied in the literature. However, the for- 
malism used here seems more convenient. We have studied the coherent and squeezed states 
of these systems in some detail. In terms of these concepts the time evolution of some wave 
packets can be discussed very conveniently. 
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NOTE ADDED AFTER PUBLICATION 

After this paper had been published we became aware of some more previous works dealing 
with the same problem [16,17]. In Ref. [16] the energy spectrum was worked out by an algebraic 
approach. In Ref. [17] the propagator was calculated by a different method from that of Refs. 
[5,6]. We thank Prof. J. M. Cervero for pointing out these papers to us. 
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